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Abstract 

The problem of reconciliation of mass balance in chemical- 
engineering systems using measured values of flows with 
incomplete data is formalized as the minimization of 
degenerate quadratic functions, with equations of mass 
balance of junctions (objects) in technological schemes. 

In order to diminish the uncertainty of data caused by 
inaccuracy of measurements and to compensate for data 
incompleteness, estimations of product flow losses done by 
experts can be introduced in the balance equations, as well 
as equations reflecting the specifics of technology, for 
example well-known relationships between feedstock and 
product flows. 

To overcome the potential degeneration and incompatibility 
in combined equations, finding normal pseudo-solutions 
was proposed for the reconciliation of mass balance. Related 
calculation formulas were explained, using matrix pseudo- 
inversion operations. 

Analysis methods of uniqueness of definition of consistent 
data were shown, as well as related absence of measure- 
ments, and detection of inconsistent constraints. A case 
study was provided. 

Keywords 
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Introduction 

The current state of theory and methods of data 
reconciliation of mass balance for industrial chemical- 
technological enterprises can be found in studies 
[Romagnoli (1999), Somov (2009), Bagajewicz (2010)] 
and in review article [Ozyrt (2004)]. Well-known up- 
to-date software applications for mass balance 
reconciliation are shown below: 

The basis of these software applications is the 
procedure of data reconciliation which can provide 
industrial enterprises with reliable information. To 
minimize reconciliation errors calculated by quadratic 
functions, least-squares method algorithms are used. 


The biggest difficulties with technological parameters 
measured values reconciliation are related to data 
incompleteness and measurement errors. 


IOO (Interactive 
On-Line Opt.) 

Academic 

Louisiana State University (USA) 
http://www.mpri.lsu.edu 

DATACON 

Commercial 

Simulation Sciences (USA) 
http://www.simsci.com 

SIGMAFINE 

Commercial 

OSI Software inc. (USA) 
http://www.osisoft.com 

VALI 

Commercial 

Belsim (Belgium) 
http://www.belsim.com 

ADVISOR 

Commercial 

Aspentech (USA) 
http://www.aspentech.com 

RECONCILER 

Commercial 

Resolution Integration Solutions 
(USA) 

http://www.ris-resolution.com 

PRODUCTION 

BALANCE 

Commercial 

Honeywell (USA) 
http://hpsweb.honeywell.com 

RECON 

Commercial 

Chemplant Technology 
http://www.chemplant.cz 


There are also: 


1) Errors in the preparation of initial data (for example, 
errors in the conversion of volume-to-mass in flows. 

2) Measurement errors or human errors (for example, 
incorrect interpretation of flow graphs of the chemical- 
technological system. 

Both types of errors can be temporarily eliminated 
during result verification. However, these errors can 
periodically reappear, for example, in the case of a 
change in the chemical composition of the feedstock or 
in the assortment of output production (a change of 
parameters in the units or in the flow graphs topology). 

Setup of the Mass Balance Reconciliation 
Tasks 

When the mass balance reconciliation takes place 
under conditions of incomplete or incorrect data, the 
following tasks appear: 

1) Actual task of mass balance reconciliation 
(getting new values for the flow rates. 
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corresponding to the mass balance equations 
and preferably close to the measured values). 

2) Diagnostics of the incomplete or incorrect data. 


The first task is a constrained minimization and can be 
examined as follows: 
f(x var ) = 



A x x var + A 2 x fix = 0 (2) 

where y = (y ; ) £ R Uvar ,z — (z ; ) £ R Uvar , are measured 
values, (y is measured at the output of some units 
(objects) or at the input of pipes, z is measured at the 
input of some units or at the output of pipes, k y — 
(kf) £ Z Uvar ,k z = (kf) £ Z Uvar , are values consisting 
of numbers 1 (if a measurement exists) or 0 (if no 
measurement exists), s y — (sf) £ R flvar , e z = (sf) e 
R n var ' are values related to measurement errors, 
x var £ R riuar , are calculated values of mass flow rates 
(variable values), x flx £ R n f ix are fixed values of mass 
flow rates (flows of feedstock, output products and 
flows measured from fluid level changes in the tanks), 
A = [A 1 |i4 2 ] = ( a t j ) £ R mxn ,n = n var + rif ix m < n , is 
the matrix of balanced equations (a,y = 1 — iq if the 
flow j enters the technological object i (unit, tank, etc.), 
a t j = -1 if flow j comes out of object i, a £ y = 0 if flow j 
does not take part in the mass balance of object i), 
u — (Uj) £ R m is the value of coefficients of loss {u L is 
the part of loss of object i from the sum of the mass 
flow rates of the incoming flows of the given object). 
Tasks (1) and (2) can be examined in conditions of 
incomplete data, that is, when the separate values y ■ 

and z • are unknown (if no measurement exists). The 
accounting of incomplete data in tasks (1) and (2) is 
done with zero k y , k z . For example, if flow rate j is 
not known, kf = kf = 0 . Incomplete initial data may 

cause negative results, such as the distortion of mass 
flow reconciliation values and the non-uniqueness of 
the solutions to tasks (1) and (2). 

Unlike incomplete data, incorrect data do not have 
clear external indicators. Incorrect initial data may 
cause negative results, such as the distortion of mass 
flow reconciliation values or the incompatibility of 
mass balance equation (2). 

Diagnostics of Incomplete or Incorrect Data 
and Mass Balance Reconciliation Methods 

Creation of an Analytical Solution to Tasks 1 and 2 
If in all measurable values of mass flow rates at least 


one of the values x j or y ■ is known (that is, all 

measurements exist), tasks (1) and (2) turn into the 
minimization of degenerate quadratic functions with 
equality constraints. The solution can be reached with 
analytical methods, for example with the method of 
Lagrangian coefficients [Crowe (1983), Romagnoli 
(1999), Narasimahan (2000), Ozyrt (2004)]. For 
example, using values c — (cy) £ R nvar , w = (wy) £ 
R n var anc J ft — (ft.) £ pn var , j-gg so l u tion X Var of tasks (1) 


and (2) can expressed by: 

x var — C -1 (w + A x (A x C~ 1 Af')~ 1 ') ■ (b — A 1 C“ 1 w)) (3) 

where 

(yy eJ ) 2 + ( ZjEf ) 2 if kj = 1 and kf = 1 
Cy=- (y jE y ) 2 if kj — 1 and kf = 0 (4) 

(zjsf) 1 if kf = 0 and kf — 1 
C = diag(c ) (5) 

yf 1 (£ ; y ) 2 4- z” 1 (fij 2 ) 2 if kf — 1 and kf = 1 
Wy = • y [ “ 1 (£ y ) 2 if kf — 1 and kf — 0 (6) 

z yl (ef) 2 if kf — 0 and kf — 1 
ft = -A 2 xf ix (7) 


The given solution does not exist in the case of 
incompatibility of balance equations (2). In addition, if 
the initial data is incomplete, formula (3) which 
contains operations of matrix inversion, cannot be 
used because the solution to the tasks (1), (2) is not 
unique. Thus, tasks (1) and (2) become the 
minimization of degenerate quadratic functions with 
equality constraints. Although solutions to similar 
tasks can be found with analytical methods, the 
solution algorithms turn out to be more complicated 
than in the degenerate cases. The following two 
solutions are possible: 

The first solution is well-known (see [1-4]) and can 
only be used when tasks (1), (2) have a unique solution, 
regardless of incomplete initial data. In this case, the 
calculations are done in three steps: 

1) The j equations containing an absence of value 
xf ar from x var are excluded from the equations 
(2), using equivalent conversions; 

2) The value x var , applicable to actual measured 
values, is calculated with formulas (3)-(7); 

3) The value x var , applicable to the absence of 
measured values, is defined from equation (2) 
taking into account the value x var that was 
calculated in the previous step. 

The second solution, presented below, consists of a 
generalized solution to tasks (1), (2) using a set of 
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matrix algebra. It does not have the defects of the 
Lagrangian coefficients and can be used for the 
reconciliation with incomplete or incorrect data, for 
example, in the case of non-uniqueness of solution of 
tasks (1), (2) caused by the absence of measurements. 


Calculations can be done with the following formula 
[Lawson (1995)]: 

x var * = Afb+(CQ - dtd 1 )) + ■ (Cw - CAfb ) (8) 

where x var is a value with n var elements(among 
which, x var , relative to the absent ones xj ar ), A* is the 
matrix, pseudo-inverse to the matrix A 1 (see [5, 6]), I is 
the identity matrix with matrix order n 1 and C, b are 
calculated with formulas (5), (7). Instead of (4) and (6), 
we now have: 


J(y/zf ) + (zy£?) if kf = 1 and kf = 1 
(yj sf ) 1 if kf — 1 and kf — 0 
( Zjsf ) 1 if kf — 0 and kf = 1 


^ 0 if kf = 0 and kf = 0 

VfHe/) 2 + z r 1 ( £ /) — land kf 

yfHe?) 2 if kf = 1 and kf = 0 

Wj = 11 

z i 1 (e ; ? ) if kf — 0 and kf — 1 
k 0 if kf = 0 and kf = 0 

Also, C — diag(c ) 


= 1 


(10) 


( 11 ) 


where c = Cy is a value with a dimension of n var , 


cj = \ Cj 1 
tO ot 


( 12 ) 


if Cj ± 0 
10 otherwise 

The validity of formulas (8)-(12) can be verified by 
examining the following task and theorem 1: 

r ) = ||Cx i;ar — Cw|| 2 -> min 
A ± x var = b 

where the symbol || - 1| is a Euclidean vector norm. 


F(x r 


Theorem 1. Tasks (1), (2) and (13) can be soluble or 
insoluble in the same sets of circumstances. If tasks (1), 
(2) are soluble, their unique normal solution (solution 
with minimal Euclidean vector norm) coincides with 
solution x var of task (13). Otherwise, the value x var is 
the unique normal pseudo-solution of tasks (1), (2) and 
(13), that is, the value with minimal Euclidean norm 
on which the minimal function f(x var ), F(x var ) is 
reached, under conditions of minimal Euclidean norm 
residuals of the equations (2). 


Calculations with Formulas (3) and (8) 

Formulas (3) and (8) contain operations of inversion 
and pseudo-inversion of matrices. 


the inversion (pseudo-inversion) of diagonal matrices 
C and C is simple and does not need any comments. 

The expression p = ■ ( b — A^^w) being 

a part of the calculation with formula (3), can be 
examined as a task of solving linear algebraic 
equations 

(A 1 C _1 A[)p — {b — Aj^-^w) (14) 

relative to the value p . Because of the structural 
specifics of matrix A 1 and also of formulas (4)-(5), the 
matrix A 1 C _1 A[ is symmetrical and positive defined. 
This condition enables the creation of iterative 
algorithms for the solution of equation (14). For 
example, it is possible to use the method of iterative 
modification of conjugate gradients which theoretically 
(in the absence of rounding-off errors) gives the 
solution of equation (14) in n var steps. 

The pseudo-inversion of matrices in formula (8) can 
also be solved by the least-squares method of some 
linear algebraic equations. However, in the case of 
incomplete initial data, the given equation can turn 
out to be degenerate, which can complicate its solution 
with well-know methods. 

Overcoming these difficulties helps to directly 
calculate the pseudo-inverse matrices, which can be 
done using the row modification of the Greville 
algorithm, as shown below [Gantmacher(1979)]. 

Greville Algorithm 

Given: matrix A of size mxn . Generally ?n A n , 
rank A < min{m,n}. 

Find: matrix A + of size n x m. 

I T 

a l a l > £ 

0 £ R n otherwise 

(cp is the first row of matrix A, Af is the first column of 
matrix A + ) 

Step k (k=2,3,...,m). 

Calculate d = a^f^. (a k is the k - th row of matrix A, 
A k _ ± is the first k - 1 of the columns of matrix A + ) 
Calculate c = dA k _ x . ( A k _ j is the first k-1 rows of matrix A) 

{ if CC T > £ 

Calculate q — l f c T 

i A k-i d otherwise 
V i +dcl T 

Calculate Af = [Aj[_ 1 — qd\ q]. 

Additional Constraints in the Task of Reconciliation 
of Mass Balance as a Way of Decreasing the Un- 
certainty of Data using A Priori Information 


The complexity of these operations is different. Thus, 


Apart from balancing equation (2), the task of mass 
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balance reconciliation may include additional relations. 
The mentioned relations can specify the structure of 
flows (for example, between units) and define 
additional balance relations (for example, between 
input and output flows of feedstock tanks) or 
proportions between feedstock input and product 
output, relative to the specific technology. 


For example, let us assume that a given feedstock unit 
(with mass flow rate x 1 ) generates a product (with 
mass flow rate x 2 ). In addition, the unit produces m 
by-products (with mass flow rates x 2 ,x 4 , ...,x m+2 ). The 
feedstock of the technological process is separated into 
main product and by-products, according to specified 
proportions that can be formalized with equations 
such as: 



*771+2 
- X ± 


0"m+2 


(15) 


Coefficients a k of the mentioned systems depend on 
the chemical composition of the feedstock and on the 
technological process. These coefficients, as well as the 
coefficients of losses iq which are present in the setup 
of tasks (1), (2), are defined by experts or calculated 
with technological process models. 


In spite of the difference of interpretation of the 
contents, the mentioned additional relations can be 
written as additional linear equations, such as: 

B 1 x var + B 2 xf ix = 0 (16) 


The integrated system of equations, composed of 
equations (2) and (16), can again be formally written as 
(2) if matrices A 1 and A 2 are redefined, and thus 

putting: A 1 


f^il 

An ^ — 

^2 

kJ 

> 

[b 2 \ 


Note that in general, the unified integrated system can 
contain tautologies (when some equations can be 
obtained from linear combinations of other equations) 
or can turn out to be incompatible (for example, 
because of a contradiction between balance equations 
and equations reflecting technological proportions 
between feedstock and product flows). 


Under the given reasons, the task of mass balance 
reconciliation with additional constraints is easier to 
solve with more universal algorithms (not sensitive to 
tautologies and incompatibility of constraints), as 
described in formulas (8)-(12), and with an analysis of 
correctness (incorrectness) of the results obtained. 


Interpretation of the Solution to Tasks (1), (2) 

Let us assume that x var is the solution to some tasks 


of mass balance reconciliation (with incomplete or 
incorrect data), obtained with formula (8). The 
following situations are examined: 

1) Let us assume xf ar , and | xj ar — yj \ » ej , 
| xj ar —Zj\»ef (if there is a discrepancy 
between calculated and measured values, 
because of inadequate measurements). 

2) Not all equations (2) have a zero difference 
between measured and real values (incorrect 
measurements). 

3) Not all values xj ar relative to the absence of 
measurements, are defined explicitly 
(incomplete initial data). 

In practice, any two or all three of the abovementioned 
situations can appear at the same time. 

The analysis of initial data to find out inaccurate 
values in cases 1 and 2 represents a nontrivial task 
which can be the subject of a separate study. 

However, in case 2 it is possible to make some express 
analyses of data, converging into the creation of a list 
of flows with fixed flow rates, linked to incomplete 
mass balance equations. This list needs to be verified. 

Analysis of Incorrect Results of Mass Balance Re- 
conciliation in the Case of Possible Incomplete 
Initial Data and Contradictory Constraints 

The incompleteness of data is not an obstacle to 
solving tasks of mass balance reconciliation. Formula 
(8) guarantees the uniqueness of all calculated values 
xj ar , for which there is at least one measured value. 
Meanwhile, values xj ar , relative to the absence of 
measurements, can still be defined uniquely in some 
cases, but not in other cases. 

To study the non-uniqueness of solutions to tasks (1), 
(2) caused by incomplete initial data, let us do the 
following: 

Assume that a matrix A 1 and a value x var split up into 
blocks: 

A 1 = [G\D],x var = [*l] (17) 

where value x G consists of those elements xj ar of the 
value x var , for which kj A 0 (that is, from measured 
values), value x D consists of those elements xf ar of the 
value x var , for which kj = 0 (that is, from non- 
measured values). Formula (17) enables you to rewrite 
equation (2) as follows: 

A 1 x var = b <=> Gx g + Dx d = b <=>' Dx D = b (18) 
where b = b — Gx g 
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FIG. 1 TECHNOLOGICAL SCHEME 


TABLE 1 MEASURED AND CORRECTED VALUES OF MASS FLOW RATES 


Flow Index 

Flow Name 

Value 

From 

To 

k y 

y , t 

s y , % 

k z 

Z ,T 

e‘ , % 

1 

Hydrofined Light Diesel 
Fraction (200-320 °C) 

*7 


Parex Process (Unitl) 

1 

1157.495 

0.8 

1 

1171.286 

0.8 

2 

Hydrogenous Gas 



Parex Process(Unitl) 

0 

- 

- 

1 

2.566 

2.0 

0 

3 

Hydrogenous Gas 


Parex Process 
(Unitl) 


1 

4.646 

2.00 

0 

- 

- 

4 

Liquid Paraffin Mixtures 
(Ranging C 9 -C 20 ) 

*7 

Parex Process 
(Unitl) 


1 

192.244 

1.08 

0 

- 

- 

5 

Dewaxed 200-320 2 C Fraction 

x7 

Parex Process 
(Unitl) 


1 

945.934 

1.14 

0 

- 

- 

6 

Hydrofined Light Diesel 
Fraction (200-320 °C) 

x7 


Parex Process (Unit2) 

1 

1341.432 

0.8 

1 

1354.663 

0.8 

7 

Hydrogenous Gas 

xT 


Parex Process(Unit2) 

0 

- 

- 

1 

8.318 

2.0 

0 

8 

Hydrogenous Gas 

x7 

Parex Process 
(Unit2) 


1 

4.660 

2.00 

0 

- 

- 

9 

Liquid Paraffin Mixtures 
(Ranging C 9 -C 20 ) 

x; m 

Parex Process 
(Unit2) 

Mixer 

1 

261.843 

1.08 

0 

- 

- 

10 

Dewaxed 200-320 2 C Fraction 

*7 

Parex Process 
(Unit2) 


1 

1065.720 

1.42 

0 

- 

- 

11 

Liquid Paraffin Mixtures 
(Ranging C 9 -C 20 ) 

*,7 

Mixer 

Charging Stock Tank 
Farm 

0 

- 

- 

1 

439.207 

5.0 

0 

12 

Liquid Paraffin Mixtures 
(Ranging C 9 -C 20 ) 

*7 

Charging Stock 
Tank Farm 

Paraffin Separation 
Unit 

0 

- 

- 

0 

- 

- 

13 

Paraffins Cn-Ci? 

*7 

Paraffin 
Separation Unit 


1 

180.631 

1.06 

0 

- 

- 

14 

Paraffins C 9 

*7 

Paraffin 
Separation Unit 


1 

13.623 

2.16 

0 

- 

- 

15 

Paraffins C 18 -C 23 

*7 

Paraffin 
Separation Unit 


1 

37.038 

2.06 

0 

- 

- 

16 

Paraffins C13 


Paraffin 
Separation Unit 


1 

41.373 

2.14 

0 

- 

- 

17 

Paraffins C 10 -C 13 

*7 

Paraffin 
Separation Unit 


1 

143.810 

2.49 

0 

- 

- 


As known from linear algebra [Albert (1972)], in the 
case of compatible equations (8), their solution can be 
described with the formula: 

x D = D + b + PAx (19) 

where 

P = I — D + D (20) 


I is the identity matrix. Ax is an arbitrary value. From 
(19) and (20) it follows that element xf is defined 
uniquely when and only when the row (or column) j 
of matrix P is entirely composed of zero elements. The 
equivalence of the given verification of rows or 
columns is caused by the symmetry of matrix P. 
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TABLE 2. FIXED VALUE OF MASS FLOW RATE 


Flow Index 

Flow Name 

Value 

From 

To 

x r " , tons 

18 

Charging Stock Tank Farm 


Charging Stock Tank Farm 


6 


TABLE 3. .LOSSES IN THE UNITS OF THE TECYNOLOGICAL SCHEME 


Unit Index 

Unit Name 

u , % 

1 

Parex Process (Unitl) 

0.35 

2 

Parex Process (Unit2) 

0.34 

3 

Mixer 

0 

4 

Charging Stock Tank Farm 

0 

5 

Paraffin Separation Unit 

0.36 


TABLE 4 INITIAL FND RECONCILED MASS BALANCES 


Component 

Before Reconciliation 

After Reconciliation 

y 

Z 

Parex Process (Unitl) 

Input, tons 

Hydrofined Light Diesel Fraction (200-320 °C) 

1157.495 

1171.286 

1157.704 

Hydrogenous Gas 

2.566 1 

2.566 

2.566 

Total: 

1160.061 2 

1173.852 

1160.269 

Output, tons 

Hydrogenous Gas 

4.646 

4.646 

4.647 

Liquid Paraffin Mixtures (Ranging C 9 -C 20 ) 

192.244 

192.244 1 

187.859 

Dewaxed 200-320 °C Fraction 

945.934 

945.934 1 

963.702 

Total: 

1142.824 

1142.824 2 

1156.208 

Loss, tons (%) 

17.237 2 (1.486 2 ) 

31.028 2 (2.643 2 ) 

4.061 (0.35 ) 

Parex Process (Unit2) 

Input, tons 

Hydrofined Light Diesel Fraction (200-320 °C) 

1341.432 

1354.663 

1342.317 

Hydrogenous Gas 

8.318 1 

8.318 

8.315 

Total: 

1349.750 2 

1362.981 

1350.632 

Output, tons 

Hydrogenous Gas 

4.660 

4.660 1 

4.661 

Liquid Paraffin Mixtures (Ranging C 9 -C 20 ) 

261.843 

261. 843 1 

253.269 

Dewaxed 200-320 °C Fraction 

1065.720 

1065.720 1 

1088.110 

Total: 

1332.223 

1332.223 2 

1346.040 

Loss, tons (%) 

17.527 2 (1.299 2 ) 

30.758 2 (2.257 2 ) 

4.592 (0.34 ) 

Charging Stock Tank Farm 

Input, tons 

Liquid Paraffin Mixtures (Ranging C 9 -C 20 ) 

439.207 1 

439.207 

441.128 

Total: 

439.207 2 

439.207 

441.128 

Output, tons 

Liquid Paraffin Mixtures (Ranging C 9 -C 20 ) 

433.207 3 433.207 3 

435.128 

Charging Stock Tank Farm 

6“ 

6 

Total: 

439.207 2 

439.207 2 

441.128 

Loss, tons (%) 

0 2 (0 2 ) 

0 2 (0 2 ) 

0(0) 

Paraffin Separation Unit 

Input, tons 

Liquid Paraffin Mixtures (Ranging C 9 -C 20 ) 

439.207 2 

439.207 2 

435.128 

Total: 

439.207 2 

439.207 2 

435.128 

Output, tons 

Paraffins C 14 -C 17 

180.631 

180.631 1 

186.096 

Paraffins C 9 

13.623 

13.623 1 

13.752 

Paraffins C 18 -C 23 

37.038 

37.038 1 

37.906 

Paraffins C 13 

41.373 

41.373 1 

41.985 

Paraffins C 10 -C 13 

143.810 

143.810 1 

158.823 

Total: 

416.475 

416.475 

433.562 

Loss, tons (%) 

16.732 (3.862) 2 

16.732 (3.862) 2 

1.566 (0.36) 


1 Value measured at the other end of the pipe (indirect measurement method) 
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2 Calculated on indirect measurement values 

3 Not measured value, but calculated from mass balance of charging stock tank farm 

4 Fixed value 


Thus, the method of verification of the non-uniqueness 
of value xf is the verification that there is a zero on 
row (or column) j of matrix P. 

In separate cases, equation (2) can turn out to be 
incompatible in case of non-reconciliation of the 
additional balance constraints and the proportions 
between feedstock and product flows, clearly specified 
in the form of additional constraints. In the given 
situation calculated by formula (8), the value x var is 
the normal pseudo-solution of tasks (1), (2): in it we 
reach the minimal function f(x var ) under conditions 
of minimal Euclidean norms errors of the equation (2). 
The equations for which the errors are not zero turn 
out to be the cause of incompatibility of equation (2). 

Therefore, to correct the incompatibility of equation (2) 
it is necessary to either exclude additional constraints 
(in relation to basic balance equations), or to correct 
them. 

Case Study 

For our case study, we will examine the results of the 
reconciliation of a part of a technological scheme in an 
oil refinery [Bannov (2000)]. See figure 1. 

The initial data for the solution of the task of 
reconciliation of the mass balance (numerical values of 
k y , k z , y, z, E y ,£ z ,u and correspondences of their 
elements to the names of flows and units of the 
technological scheme) are shown on tables 1-3. The 
character in tables 1 and 2 indicates other objects (units) 
of the technological scheme of the oil refinery that are 
not shown on figure 1. The character - in tables 1 and 2 
indicates the absence of relative measured value. 

Values of initial and reconciled mass balances are 
shown on table 4. 

This example is characterized by the presence of 
uncertainties (no measurements for liquid paraffin 
mixtures) from the charging stock tank farm to the 
paraffin separation unit. The model of balance 
reconciliation contains an additional constraint, but it 
does not result in data incorrectness. 

In this example, the tool of data reconciliation is the 
generalized solution to the tasks (1), (2), based on 
formulas (8)-(12). The balance reconciliation has only 
one solution, as shown below. 

An additional constraint is the dependence between 


the mass flow rates of liquid paraffin mixtures 
appearing in the paraffin separation unit ( x™ ) and 
the mixture of paraffins C 10 -C 13 ( + x ™ ): 

0A5x^ r - x^ r - xtf r = 0 ( 21 ) 


The matrix and right part of the balance equation with 
additional equation (21) is as follows: 


( 0.9965 

0.9965 -1 - 

1 -1 0 

0 0 0 

0 0 0 

0 

0 

0 

0 

0 ) 

0 

0 0 

0 0 0.9966 

0.9966 -1 -1 

-1 0 0 

0 

0 

0 

0 

0 

0 

0 0 

1 0 0 

0 0 1 

0 -1 0 

0 

0 

0 

0 

0 

0 

0 0 

0 0 0 

0 0 0 

0 1 -1 

0 

0 

0 

0 

0 

0 

0 0 

0 0 0 

0 0 0 

0 0 0.9964 

-1 

-1 

-1 

-1 

-1 

0 

0 0 

0 0 0 

0 0 0 

0 0 0.45 

0 

0 

0 

-1 

- 1 , 


0 


u — 

6 

0 

A 

In accordance with formulas (8)-(12), 


' 0.1 52 3 


' 26.839' 


' 6.586' 


6 1 157.704' 

19.486 


974.279 


0.051 


2.566 

10.762 


538.097 


0.093 


4.647 

0.482 


44.596 


2.076 


187.859 

0.093 


8.134 


10.784 


963.702 

0.131 


23.182 


7.625 


1342.317 

6.011 


300.553 


0.166 


8.315 

10.730 


536.481 


0.093 


4.661 

0.354 

, w = 

32.742 

, c = 

2.828 

,x var * = 

253.269 

0.066 


4.654 


15.133 


1088.110 

0.046 


0.911 


21.960 


441.128 

0 


0 


0 


435.128 

0.522 


49.272 


1.915 


186.096 

3.398 


157.333 


0.294 


13.752 

1.311 


63.624 


0.763 


37.906 

1.129 


52.778 


0.885 


41.985 

, 0.279 ) 


v 11.215, 


, 3.581, 


[ 153.823, 


In accordance with (17), matrix D consists of one 
column: column 12 of matrix A lr corresponding to non 
-measured value r (mass flow rate from the charging 
stock tank farm to the paraffin separation unit): 

( 2 \ 

D = -1 

0.9964 
\ 0.45 ) 

Since D is a column, matrix P, created in accordance 
with (20), degenerates into a scalar, therefore P — 0, 
from which follows the uniqueness of definition of 
value x\? r ' . 
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The substitution of value x var into the equation 
A ± x = b (balance equation subject to additional 
equation (21)) shows that the error turns out to be 
equal to zero, which confirms correctness of data and 
results of the solution to the task. 
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